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ABSTRACT 

Tidal disruption events (TDE) in which a star is devoured by a massive black hole at a galac¬ 
tic center pose a challenge to our understanding of accretion processes. Within a month the 
accretion rate reaches super-Eddington levels. It then drops gradually over a time scale of 
a year to sub-Eddington regimes. The initially geometrically thick disk becomes a thin one 
and eventually an ADAF at very low accretion rates. As such, TDEs explore the whole range 
of accretion rates and configurations. A challenging question is what the corresponding light 
curves of these events are. We explore numerically the disk luminosity and the conditions 
within the inner region of the disk using a fully general relativistic slim disk model. Those 
conditions determine the magnitude of the magnetic field that engulfs the black hole and this, 
in turn, determines the Blandford-Znajek jet power. We estimate this power in two different 
ways and show that they are self-consistent. We find, as expected earlier from analytic argu¬ 
ments (IKrolik & PirarJ 120121) . that neither the disk luminosity nor the jet power follows the 
accretion rate throughout the disruption event. The disk luminosity varies only logarithmi¬ 
cally with the accretion rate at super-Eddington luminosities. The jet power follows initially 
the accretion rate but remains a constant after the transition from super- to sub- Eddington. At 
lower accretion rates at the end of the MAD phase the disk becomes thin and the jet may stop 
altogether. These new estimates of the jet power and disk luminosity that do not simply follow 
the mass fallback rate should be taken into account when searching for TDEs and analysing 
light curves of TDE candidates. Identification of some of the above mentioned transitions may 
enable us to estimate better TDE parameters. 

Key words: accretion, accretion discs - black hole physics - relativity - methods: numerical 
- galaxies: jets 


1 INTRODUCTION 


The surprising X-ray emission from the tidal disruption events 
(TDE ) candidates Swift J1644 jBloom et al] 2011: Leva n et al.l 
1201 ll:lBurrows et all201 lh and Swift .12058 dCenko et al.ll2012h led 
to a re-examination of the accreti on pro cesses that take place within 
these events. 25 years ago. lReesI jl988l) outlined the basic dynami¬ 
cal processes relevant to TDEs. A star is disrupted by a supermas- 
sive black hole. The stellar material is spread out and it returns, 
after a delay, to the vicinity of the black hole where it forms an 
accretion disk. If the disrupted stellar material has a uniform dis¬ 
tribution of orbital binding energy per unit mas£], the infall rate 
onto the accretion disk satisfies: dM/dt oc f ' 5/3 jph inne vi fl989h . At 
the peak accretion rate, the luminosity would be super-Eddington 


* E-mail: tsvi.piran@mail.huji.ac.il (TP); asadowsk@mit.edu (AS); 
atchekho@berkeley.edu (AT) 

1 This assumption that was verified by various numerical simulations, see 
e.g. Lodato. King & Pringlel f2009l) ;lGui llochon & R amirez -Ruizl j2013l) 


{yT nieill 19991) . Within this disk, gravitational energy is dissipated 
and the heat is radiated in the usual quasi-thermal fashion. If it is 
thermally radiated by the accretion disk, the associated tempera¬ 
ture would be in the extreme ultra -violet (EU V) or perhaps the soft 
X-ray band (see e.g. iLodato & Rossill201ll) . Consequently, most 
searches for such events hitherto have been carried out in the EUV 


region. Indeed several such candidates were found (Gezari et al. 

2008 

; Cappelluti et al. 20091; Ivan Velzen et al. 2011; Arcavi et al. 

2014 

). 


In early 2011, iGiannios & Metzgeij (1 201 ll) suggested that 
TDEs would involve ejection of relativistic jet that would give rise 
to a radio signal when it is slowed down by the surrounding matter. 
This prediction was verified when the radio emission from Swif t 
J1644 was discovered a few months later JZauderer et cil. imj. 
Howe ver, it turned out that the jet is also a p owerful source of X- 
rays jBloom et al,ll201 ll : iBurrows et al,ll201 lh . These X-rays arise, 
most likely, from the inner part of the jet, probably due to internal 
shocks taking place there. A few months later. ICenko et al.l d2012h 
discovered a similar X-ray signature from another TDE candidate 
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(Swift J2058). The latter source also showed an optical emission, 
most likely the thermal component arising from the accretion disk 
llCenko et al 1 2012h . 

The interplay between the (unexpected) nonthermal jet emis¬ 
sion and the thermal disk component is intriguing. Over a short 
period of a few months the accretion rate spans a large range of 
values, beginning at super-Eddington, turning to sub-Eddington 
power and then diminishing. IKrolik & Piranl 1201 21) addressed this 
issue using simplified models for both the thermal luminosity and 
for the jet power. For the former they assumed that the luminos¬ 
ity is capped at the Eddington luminosity at the super-Eddington 
phase and then it decreases following the mass accretion rate. 
For the jet power IKrolik & Piran j 20 12h related the Blandford- 
Znajek (hereafter BZ, iBlandford & Znaiekl 1 97~7l : lMcKinne\l2005l ; 
lElawlev & Krolikll 20061) jet power to the pressure at the inner parts 
of the disk. This pressure is linearly related to the mass accretion 
rate at the super-Eddington phase. But the pressure is independent 
of the mass accretion rate during the r adiation dominated phase , 
that foll ows the super-E ddington phase (Moderski & Sikorai fl996h . 
This led IKrolik & PiraJ d2012t) to suggest that the jet luminosity 
will be constant during this phase. At sufficiently low accretion 
rates the disk becomes gas pressure dominated and the pressure 
will decrease as the accretion rate decreases. Thus one would ex¬ 
pect two transitions in the jet’s and disk’s light curves. One when 
the accretion rate drops below Eddington and the other when the 
disk becomes gas pressure dominated. A third transition would 
arise at very low accretion rate, once th e disk becomes an advec- 
tion dominated accr etion flow (ADAF) dAbramowicz et alii 19881 : 
iNaravan & Yil 1 9951) . 

Our goal here is to confront these rather simple estimates 
with a more detailed computation of the accretion disk structure 
as well as with recent numerical results on the jetOuminosity. To 
this, end we use relativistic slim disk solution JSadowsk~ill20l ll) to 
estimate the accretion disk structure over a large range of accre- 
tion rates. The slim disk m odel generalises the standard thin disks 
dNovikov & Thornell 1 9731) towards high accretion rates. It allows 
for non-Keplerian rotation, radial gradients of pressure, and advec- 
tion of heat. As a result, when the accretion rate is super-Eddington, 
the disks are no longer radiatively efficient and their angular mo¬ 
mentum profile can significantly depart from Keplerian. The slim 
disk model uses a viscosity and assumes the accretion rate is con¬ 
stant, i.e., there are no outflows. 

Simulations of jet formation from accreting BHs have led 
to different estimates of the jet luminosity. The early, pioneer¬ 
ing simulations of jet formation in ge neral relativistic magnetohy- 
drodynamic_(GRMHD) simulations dDe Villiers & Hawlevl 120031: 
lHawlev & Kroliki 120061) suggested that the magnetic pressure in 
the funnel regions of the jet is determined by the thermal pres¬ 
sure of the accretion disk near the BH horizon. It also ap¬ 
peared that the accretion disk structure did not seem to be af- 
fect ed by the prese nce of the jets or by how strong the jets 
are dBeckwith et al.ll2008l) . Later, it became clear that both disk 
thermal pressure and jet magnetic pressure can be much higher 
than previously simulated in GR a nd can be as high as the 
ram pressure of the infalling g as (|Tchgkhoyskoy_etalJ 1201 ll : 
ITchekhovskov & McKinnevI 1200 : iMcKinnev et al. 2012lf . When 
the jet pressure reaches this limit, the magnetic field is strong 
enough to obstruct the accretion of gas onto the BH, and this 
leads t o the formation of a ma gnet ically arrested disk, or a 
MAD ( ITchekhovskov et ail 1201 ll : see llgumenshchev et akl 1200 3l : 
Ilgumensh chevl20 0^for sim ulations of MADs in Newtonian relativ¬ 


ity and Bisnovatvi-Koean & Ruzmaikinlll 9741 . 119761 : INaravan et al.l 


120031 : ITchekhovskov 2 0151 for an analytic consideration). In the fol¬ 
lowing we will combine the simpler jet luminosity estimates based 
on the disk pressure with the model fits coming from sophisticated 
numerical simulations. 

We begin in ^examining different methods for estimating 
the jet power. In @we discuss the slim disk model of an accretion 
flow. We present both the thermal emission and the jet power and 
we compare the expected behaviour of the jet and disk luminosity 
for different accretion rates in 0 We examine the implications for 
these findings to tidal disruption events in ^5] We conclude in ^6] 
with a discussion of the limitations of the analysis and possible 
observational implications. 


2 JET POWER 

At the order of magnitude level, the Poynting luminosity of a BZ 
jet that emerges from the vicinity of a black hole can be estimated 
through a simple dimensional argument. The local magnetic energy 
density near the black hole is proportional to B 2 , where B is the 
poloidal field intensity near the horizon. The jet power depends on 
this energy density and on the area from which the jet emerges. 
This area is of order >;rr 2 , where r g is the gravitational radius of 
the black hole. Combined the overall luminosity can be written as 
^JC, = f(a/M) c (B 2 /&7t) (ttr 2 ), where M is the black hole’s mass 
and a, = a/M its specific angular momentum. The function f(a ,) 
is dimensionless and we approximate it here as f(a,) « a 2 (e.g., 
ITchekhovskov et aHl20l 2l : lTchekhovskovll20 1 5l) . The task now is 
to better estimate the magnetic field ( B ) and the area from which 
the outflow emerges. 

2.1 The pressure formula 

The jet is in pressure balance with the inner parts of the accretion 
disk, therefore the strength of the magnetic field is determined by 
the pressure in the inner disk. iBeckwith et al .1 J2008I) have demon¬ 
strated that the magnetic pressure near the horizon is generally 
bounded above by the maximal pressure in the equatorial plane 
near the inner edge of the disk, p max , and it is bounded below by 
the magnetic pressure at that location. Thus, we can estimate the 
strength of the magnetic energy density using p max , the maximal 
pressure near the edge of the disk: 

B 2 /8n = Pmax, (1) 

where the factor f} B is an unknown dimensionless factor of order 
unity. A second factor is the size of the region from which the jet 
emerges. We use the radius of the innermost stable circular orbit 
(ISCO), Risco. t° characterise the size of the inner disk, which is 
comparable but smaller than the position of the maximal pressure. 
Together, we can write the BZ luminosity as: 

Rjci.p = /f/;^R[sco/Aw»AA ~ 7r ^’rsco/-’inax< 2 ;f'. (2) 

The unknown coefficient will be determined in Section [23] by 
comparing with the jet power obtained in numerical simulations. 


2.2 The H/R formula 


As in the previous section, we will balance the magnetic pressure 
pressure in the jet funnel against the pressure of the accretion disk. 
The vertical force balance approximately gives. 


- * 0,1 H 1 = 


R 2 


(3) 
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where p and p, are the midplane pressure and the corresponding 
density of the disk, respectively, I2 K and Vr are the Keplerian ve¬ 
locities at radius R, and H is disk half-thickness. Applying the ver¬ 
tical equilibrium at the radius of the pressure maximum, R, to eq.[2] 
we get, 

P ia ,niR=d>Vl^Halc, (4) 

where we introduced the surface density E = 2 pH, and allow the 
coefficient P' B to dilfer from fi B . The rest mass conservation re¬ 
quires, 

M = htSLV t , (5) 


where M is the accretion rate and V r is the absolute value of the 
radial velocity. Using this formula we get, 


Pb 


vl H 


r jet H/R = -Me 2 —-a;. 
jet,H/K 2 y R , 


(6) 


For radiatively inefficient accretion flows, the magnitude of the ra¬ 
dial velocity V r is comparable to the Keplerian velocity Vk and does 
not depend on the accretion rate. We may therefore approximate in 
this case eq.[6]and write, 


Fjet,H/R Me 2 —a 2 ,. 

K 


(V) 


2.3 Normalisation 


Numerical simulations of jets in bot h optically thin and thick , 
radiatively inefficient MAD disk s jTchgkhovskoy etalj 


Tchek hovskov & M cKinnevi 120121 : iMcKinnev et al 


2012 


provide the missing scaling factor and give (ITchekhovskov 



jrj _ . iH/R 2 

Fjet,H/R,RIAF - 1.3 MC «„• 


0.3 


( 8 ) 


Because of the taken assumptions, this formula for the jet power 
is valid only for radiatively inefficient disks — ADAFs or super¬ 
critical (or super-Eddington) disks, i.e., in the limit of lowest and 
highest accretion rates. To estimate the power of the jet in disks 
with not so extreme accretion rates, one has to use more general 
formulation (eq.[2]or eq. [6j. 

All the three formulae are expected to give the same estimate 
of the jet power. To satisfy this condition, we choose the coefficients 
Pb (eq. 0 and p' B (eq. 0 so that the jet power estimates agree with 
eq.[8]for M » Mem, and finally get, 


and 


Fjet.p — 4.0 nR lsco p mdx aic, 


. n Vl H/R n 
^jet.H/R = 0.16Mc 2 —^-—3—aj. 


V r c 0.3 


(9) 

( 10 ) 


3 DISK MODEL 

To model an accretion disk we use the general relativistic slim disk 
solutions of Sad owski l J20 1 1|). The slim disk model generalises the 
stand ard thin disk dShakura & Sunvae~vlll973l : [Novikov & Thome! 
1 197.31) to arbitrary accretion rates. It allows for non-Keplerian ro¬ 
tation and advective cooling by photons trapped in the flow. It as¬ 
sumes constant accretion rate, i.e., it does not allow for outflows, 
and it adopts the a prescription for viscosity. The slim disks reduce 
in the limit of small accretion rates to the standard thin, relativistic, 
Keplerian disk. 



Accretion rate [Edd.] 

Figure 1. The maximal pressure in the disk as a function of the accretion 
rate for BH spins a, = 0.3, 0.6 and 0.9 and M = 10 1 M a . 


The system has a characteristic luminosity, the Eddington lu¬ 
minosity, L E dd'- 


AnGMc M 

i-Edd = -= 1-25 x 10 38 — erg/s. 

k,. M a 


( 11 ) 


This corresponds to the Eddington accretion rate, here defined as: 


1 U 


Mr., 1,1 —- 


1 4 tcGM 


= 2.44 x 10' 


Mr. 


■ g/s, (12) 


rj c rj CK es 
where we put the efficiency of a thin disk around a non-rotating 
BH, r\ = 0.057, and /c es = 0.4cm 2 /g. Once the accretion rate is 
near and above MEdd. photons do not have enough time to diffuse 
out of the disk and a fraction of them is adverted with the flow. 
This extra advective cooling modifies the structure of the disk. In 
partic ular, the disk radiates less efficiently and its luminosity scales 
as JPaczvnskill 198(1) . 


L*L Ed d(l+log-^-). (13) 

\ Mraa I 

The power of the jet depends on the parameters of the un¬ 
derlying accretion flow. The pressure formula for the jet power 
(eq. |9j is parametrised in terms of the maximal total pressure in 
the equatorial plane, p max . Fig- 0 presents p mdx as a function of 
the accretion rate for three values of BH spin: a, = 0.0 (thickest), 
0.6, and 0.9 (thinnest line), assuming a BH mass M = 10 1 M Q . For 
the lowest accretion rates, M < O.lMEdd, the disk is gas pressure 
domi nated and the pressure at fi xed radius is expected to follow 
M 4/5 dShakura & Sunvaevll 19731) . However, the radius of the pres¬ 
sure maximum is not fixed and the profiles of this quantity follow 
this dependence only qualitatively. When accretion rate exceeds 
~ 0. lMsdd the disk is radiation pressure dominated. The standard 
thin disk theory predicts that the pressure at the equatorial plane 
at a given radius is independent of the accretion rate. This explains 
the flattening of p max profiles around 0.1A/Edd- Because the pressure 
maximum is not at a fixed radius and because the advective cool¬ 
ing gradually kicks in when approaching Mem, the p max slightly 
varies with M. For super-critical (exceeding M em ) accretion rates, 
the disk enters the slim disk regime and the maximal pressure is 
proportional to the accretion rate. 

The H/R -based estimates of the jet power (eqs. [8] and ITOt 
depend on the disk thickness which we parametrise by the maxi¬ 
mal disk opening angle, 0 h = arctan (H/R). Fig. [2] presents pro¬ 
files of disk thickness for various accretion rates. BH spins and 
masses. The top panel corresponds to a non-rotating BH. For m = 
M/Mem = 0.3 the thickest radiation pressure dominated region 
is located between R = 10 and 100. Disk thickness is determined 
there by the local radiative flux. Because both the flux and the ver¬ 
tical component of the gravity are proportional to the BH mass, the 
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Figure 2. Radial profiles of disk thickness (defined as the opening angle 
©H = atan (H/R)) for various accretion rates. The top panel compares two 
BH masses (10 6 vs 10 s M a ) assuming zero BH spin, while the bottom one 
compares two values of BH spin (a, = 0 vs 0.9) for a BH mass 10 7 M Q . 

disk opening angle does not depend on the BH mass. Further out, 
where gas pressure dominates, the BH mass has a slight impact on 
disk thickness. The extent of the radiation pressure dominated re¬ 
gion increases with accretion rate and reaches R > 1000 for highly 
super-critical accretion rates. 

The bottom panel of Fig.[2]depicts the disk opening angle for 
a fixed BH mass M = 10 7 M Q and two values of BH spin, a , = 0.0 
(green) and 0.9 (blue lines). The disks (and hence the thickness pro¬ 
files) extend more inward for the rotating BH. This reflects the fact 
that the radius of the innermost stable circular orbit decreases with 
increasing BH spin. At the same time, for a fixed accretion rate, 
disks around rotating BH have a higher maximal thickness than 
their non-spinning counterparts. This reflects the higher accretion 
efficiency that results in a higher luminosity. 

Finally, in Fig.[3]we plot the maximal disk opening angle as a 
function of accretion rate for three BH spins. When the disk is gas 
pressure dominated at all radii (M < 0.001 Me<m) the maximal disk 
thickness is located far from the BH and it does not depend on its 
spin — the lines therefore coincide. In the intermediate, radiation 
pressure dominated, regime the disk thickness increases with the 
BH spin, as discussed above. For the radiatively inefficient (slim) 
regime, the maximal disk thickness saturates at 0 H ~ 0.3. This re¬ 
flects the fact that the thickness of advection dominated accretion 
flows depends only on the ratio of the radiation to gas pressure, and 
approaches the limiting va lue 0 H ~ 0.3 for radiation pressure dom - 
inated super-critical disks dNaravan & Yil 19941 : 1 Vieira et al.ll2015l) . 


4 JET AND DISK LUMINOSITIES 

The standard model of thin disks dShakura & Sunvaevll 19731) pre¬ 
dicts that energy liberated is proportional to the mass accretion, M, 
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Figure 3. The maximal opening angle, ©h, as a function of the accretion 
rate for BH spins a, = 0.3, 0.6 and 0.9. 


and determined by the accretion efficiency, 77 , 

L = r/Mc 2 . (14) 

For non-rotating BH 77 = 0.057. As discussed above, once the ac¬ 
cretion rate approaches and exceeds the Eddington limit, the effi¬ 
ciency is decreased. These facts are reflected in the disk luminosity 
profile shown in Fig. [4] corresponding to a BH spin a, = 0.6 and a 
BH mass Mbh = 10 1 M Q , with the red line. As long as M < Me dd, 
the disk luminosity increases linearly with the accretion rate. Once 
this limit has been exceeded, the luminosity grows roughly with 
logarithm of M (eq. E3. 

We now apply the jet power formulae derived in Section[2]and 
estimate the jet power for each accretion rate using the correspond¬ 
ing slim disk solution as the underlying disk model. 

The pressure formula (eq. [9j reflects the fact that the disk and 
jet pressures balance each other. Knowing the disk thermal pressure 
we may therefore estimate the jet magnetic pressure, and there¬ 
fore the magnetic flux in the jet. This quantity, together with the 
known BH spin, provides an estimate of the jet power. Its depen¬ 
dence on mass accretion rate is shown in Fig. [4] with the blue line. 
For a fixed BH mass and spin eq. [9] depends only on the maxi¬ 
mal value of the disk thermal pressure, and therefore the profile 
of the jet power estimated this way resembles the profile of the 
corresponding maximal disk pressure (Fig. Q}. The jet power in 
the super-Eddington regime grows proportionally to the accretion 
rate, remains roughly constant in the radiation-pressure dominated 
thin disk regime ( 0.1 < M/Meaa < 1), and follows M 4/s for gas- 
pressure dominated thin disks. These properties result in very pow¬ 
erful jet power estimates (exceeding the disk luminosity by 2-3 or¬ 
ders of magnitude) for the thinnest disks. This is inconsisten t with 
observations (see, e.g., Fender et al. 2004; Russell et al.ll201 lh and 
suggests that the assumptions behind the pressure formula break 
down in this regime. 

The other formula for the jet power (the H/R formula, eq. ED, 
similarly uses the disk thermal pressure as the proxy for the jet 
magnetic pressure and magnetic flux at the horizon. However, the 
disk pressure is replaced with the disk thickness using the vertical 
equilibrium equation, and the radial velocity V r and the disk thick¬ 
ness H are introduced. These are taken directly from the numerical 
solutions of slim disks. The corresponding jet power is plotted with 
the green line in Fig. [4] It coincides with the pressure formula for 
M > Meaa and stays close for lower accretion rates. This fact proves 
that both formulae properly identify the disk thermal pressure, al¬ 
though the H/R formula does it indirectly. From now on we will 
not distinguish between these two ways of estimating the jet power 
and use eq.^as the proxy. 
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Figure 4. Jet (green and blue) and disk (red line) luminosities for Mbh = 
10 1 Mq and BH spin a* = 0.6 a function of accretion rate. The green line 
corresponds to eq. urn The blue line shows the jet power calculated accord¬ 
ing to eq.[9] which we take below as our fundamental jet power estimate. 


In Fig. [5] we plot the jet power and its efficiency, 


•'-'jet 

m?-’ 


as a function of accretion rate for three values of BH spin with 
solid blue and dotted black lines, respectively. The profiles of the 
jet power for all spins agrees with the discussion above. The jet 
power is expected to grow roughly with BH spin squared, there¬ 
fore, the higher the spin, the larger is the jet power. In the limit 
of highest accretion rates, the efficiency calculated according to 
both formulae for spin a, = 0.9 is close to 77 = 1, what re¬ 
flects the fact that the efficiencies of jets observed in GRMHD 
simulations may be very high, and even exceed l_for the highest 
spins and magnetically arrested d isks dTchekhovskov et akll201ll : 
iTchekhovskov & McKinneJl2012h . 

The jet efficiency remains roughly constant as long as we are 
in the radiatively inefficient regime ( M > Meaa). Below that limit 
the efficiency increases with decreasing accretion rate. Because 
the disk thermal pressure stays roughly constant in the radiation- 
pressure dominated, thin disk regime, when the accretion rate 
changes by an order of magnitude, or more, the efficiency must 
increase by that amount. Furthermore, the pressure dependence on 
accretion rate for gas-pressure dominated disks is P oc M 4/5 . what 
further increases the efficiency. Such extreme efficiencies for the 
thinnest disks (rj ~ 100 for 0 . 001 M Edd ) are unphysical and incon¬ 
sistent with the lack of observed jets in high/soft state of Galactic 
X-ray binaries. 

Which assumptions we took break down in the sub-Eddington 
regime? Why are the jets not there? It is probably the assump¬ 
tion that p B that determines the ratio of the magnetic pressure to 
the disk pressure remains constant and is independent of the dra¬ 
matic changes that take place in the disk when it transits from 
a thick to thin disk and from radiation pressure to gas pressure 
dominated regime. In particular, one expects the radial disk ve¬ 
locity to decrease and the large-scale magnetic flux to diffuse 
outwa rd faster than it is being advected inward by the accretion 
flow iLubow et al.lll994l: iGuilet & OgilviejEoIl. 120131: see how¬ 
ever Rothstein & Lovelace 2008). In addition, one expects for thin 
disks a larger gap between the inner part of the disk and the black 
hole. As the accretion rate decreases, less and less material is within 
this gap and the thin disk may not be able to drag the necessary 
magnetic field to the vicinity of the black hole. 
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Figure 5. Jet power (blue lines, left y— axis) and its efficiency (given by 
eq. 0 and shown with dotted black lines, right y-axis) for Mbh = 1O 7 M 0 
and BH spins a , = 0.3, 0.6 and 0.9. 


5 JET AND DISK IN TIDAL DISRUPTION EVENTS. 

We turn now to apply the previous result to a TDE. We consider a 
star with mass M, (measured in solar masses) that is disrupted by a 
massive black hole (BH). We approximate the main sequence mass- 
radius relation by R, ~ £ =; 0.2_for 0.1 < M, < 1, but 

increases to ^ 0.4 for 1 < M, < 10 {iKippenhahn & Weigertll 19941) . 
Finally we define k as the apsidal motion constant (determined by 
the star’s radial density profile) and / is its binding energy in units 
of GM;/R, {Ph innevl 19891) . In numerical estimates, we scale k/f to 
the value for fully radiative stars, 0.02, because this is a reasonable 
approximation for main sequence stars with O.4M 0 < M < 10M e 
dKippenhahn & Weigerll994l) . 

Assuming that the mass accretion rate follows the fall back 
time of stellar material onto the central black hole, and that the 
disrupted star has a uniform distribution in orbital binding energy 
per unit mass, matter returns to the region near the pericenter ra¬ 
dius at a rate M oc (r/f 0 )~ 5/3 lPhinnevl ll989i) . The characteristic 
timescale to for initiation of this power-law accretion rate is_the 
orbital period for the most bound matter dLodato, Kin g & Pringle 
120091: ItCrolik & Piranll2012h : 

to “ 15 x 10 5 itff- 3 « /2 M^ 7 s. (15) 

Using this time scale we calibrate the maximal accretion rat^jas: 

44peak ~ 0.3 x 10 27 Mi 1+3f,/2 M g ^(fe//)- 1/2 g s" 1 . (16) 

We convolve now the previous estimates of jet luminosity and 
disk power with the accretion rate evolution to obtain the expected 
light curves. These are shown for 10 6 M e and 10 1 M e BHs with 
a, = 0.6 spin in Fig. [6] A quick inspection reveals that, as ex¬ 
pected, even with the more detailed calculations, the thermal disk 
lig ht curves (red l ines) follow more or less the simple estimates 
of llCrolik & Piranl d2012t) . The disk luminosity is roughly constant 
at short time scales (when M > Mem). After the transition to 
M < Mem , which takes place roughly at t/t 0 = 6 and 20 for 10 1 M e 
and 10 6 M o BHs, respectively, the disk luminosity decreases pro¬ 
portionally to the mass accretion rate. For a higher BH mass (dot¬ 
ted lines) the qualitative behaviour remains the same. However, be¬ 
cause the initial accretion rate (eq. EH is lower, the disk enters ear¬ 
lier the sub-Eddington and thin disk phases. 


2 See however Ishiokawa et IDIml) for caution concerning the onset of 
accretion in TDEs and the possibility of a lower maximal accretion rate that 
takes place at a later moment. 
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t/to 


Figure 6. Disk luminosity (red) and jet power (blue line) as a function of 
time for a tidal disruption event of a Af* = 1 star by a Mbh = 10 6 M Q (solid) 
or 10 7 M o (dotted lines) BH. Other parameters: M* = 1, k/f = 0.02, and 
£ = 0 . 2 . 

The jet lumino sity ev olution is also consistent with predictions 
of iKrolik & Piranl d2012h . Initially, as long as the disk is super- 
Eddington, the jet power (denoted with blue lines) decreases pro¬ 
portionally to the accretion rate. Once it enters the sub-critical, ra- 
diatively efficient regime, the jet power hardly changes reflecting 
the fact that the disk pressure is not sensitive to the accretion rate. 
This changes once the disk becomes gas pressure dominated. Our 
simplistic formulae predict further decrease of the jet power, al¬ 
though not as steep as in the super-critical regime. However, as dis¬ 
cussed earlier, thin disks are not likely to sustain strong magnetic 
flux at the BH and therefore should produce less efficient jets than 
we predict. 


6 DISCUSSION AND CONCLUSIONS 

We have used numerical modelling to estimate the jet power and the 
accretion disk luminosity as a function of an accretion rate over a 
very wide range of accretion rates. Our accretion disk model use the 
relativistic slim disk solutions of lSadowskil d201 ll ). It switches from 
a thick disk for accretion rates above Eddington to a thin disk be¬ 
low. We have used two different expressions for the jet power, one 
based on the maximal pressure and another, motivated by numeri¬ 
cal simulation, based on the H/R ratio. Remarkably, both formulae 
give comparable results. 

Our results confirm earlier expectations dKrolik & Piranl20l2h 
that were based on simple scaling arguments. According to these 
expectations the jet power and the disk luminosity do not simply 
follow the accretion rate. Specifically, for super-Eddington accre¬ 
tion the jet power follows the accretion rate but the disk bolometric 
luminosity varies only logarithmically with it. At sub-Eddington 
accretion the jet power becomes roughly constant while the disk 
bolometric luminosity follows the accretion rate. A puzzling be¬ 
haviour in our solution is the appearance of a very large efficiency 
(jet power vs. accretion rate) at low accretion luminosity. While 
puzzling we note that this is not impossible in principle as the 
BZ jets gets its power from the rotational energy of the central 
black hole and as such could, in principle at least, operate even 
with lower accretion rate (provided that it is surrounded by strong 
enough magnetic fields). However, as we discuss shortly below, it 
is likely that the jet power switches off completely in this low ac¬ 
cretion rate regime. 

It is interesting to note that in our model the jet power is higher 
by at least one order or magnitude than the disk luminosity. This 
holds as long as a jet exists for the whole range of accretion rates 
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Figure 7. Similar to Fig.[S]but only for Mbh = 10 6 M Q . The labels above 
the top axis denote the accretion rate in Eddington units at given time. The 
plateau before t = 3fo schematically reflects the pre-MAD stage when the 
accretion rate is too large to form magnetically arrested disk. The jet power 
was arbitrarily damped for t > 200fo to reflect the fact that the disk enters 
the thin disk regime (M < 0. IMEdd) and is not expected to produce efficient 
jet at all. The dotted blue line reflects corresponding decrease in the jet 
power if the transition takes place already at M^dd . 

and time scales considered. This may not be related directly to the 
question whether it is easier or more difficult to detect the jet emis¬ 
sion or the disk emission as those depend also on the spectral en¬ 
ergy distribution of these two components and on sensitivity and 
coverage of different detectors in the relevant spectral regimes. 

It should be stressed that three possibly important effects have 
not been taken into consideration in our work. First, we assumed 
(implicitly) that the magnetic fields needed for the BZ mecha¬ 
nism to operate are there right from the beginning. This might 
not be the case in TDEs in which the infalling material is not 
highly magnetised initially. 0 Thus it may take some time to build 
this magnetic field and in this case w e expect a constant (or eve n 
increasing) jet power at early times dTchekhovskov et aI1l2014l) . 
This effect was marked in Fig. [7] as a solid horizontal line marked 
“pre-MAD”. In other transient astrophysical source, core-collapse 
gamma-ray bursts, this stage lasts for the most duration of the event: 
in fact, then the MAD_onset plausibly marks the end of the burst 
dTchekhovskov & Gianniosll2015l) . 

A second potential effect that we have ignored here is the pos¬ 
sibility of the jet switching off when the disk makes a transitio n 
from thick to thin dZauderer et al.l2013l : lTchekhovskov et al.l2014l) . 
This transition takes place shortly after the accretion rate makes the 
transition from super to sub-Eddington. As mentioned earlier such 
a turn off of the jet can arise from the fact that the weaker accretion 
rate does not push the magnetic field sufficiently close to the black 
hole, or that the large-scale magnetic flux diffuses outward more 
effectively than it is dragged inward by the disk. This possible ef¬ 
fect was marked in Fig. [7] by arbitrary damping the jet power for 
M < 0.IMEdd what indicates that the jet power may not satisfy the 
simple curve in this region. Note that in the future, as mass accre¬ 
tion rate drops even further, M < 0.01M E dd, the accretion disk is ex¬ 
pected to transition to a geometrically-thick radiatively-inefficient 
ADAF, which can c ause re-launching of the j ets and X-ray emis¬ 
sion from Sw 1644 dTchekhovskov et alj2014l) . 

Finally, we assumed that all of the mass fallback ends up 

3 The energetics of Sw J1644 requires the presence of a large-scale mag¬ 
netic flux that exceeds by 3—5 orders of magnitude the magnetic flux 
through a main-sequence star. Such a large magnetic fl ux could be captured 
by the debris stream from a pre-existing accretion disk iTch ekhovsko y et alj 
|2014|) , as recently seen in numerical simulations jKelley et al i2014l) . 
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reaching the black hole. The processes of gas circularisation are 
not well-understood, and it is possible that the fraction of gas that 
reaches the hole depends on the Eddington ratio, and this can cause 
additional deviations from simple power-law scalings in time. Note 
also that the scaling vs time of intensity in a particular detector 
might deviate from a power-law due to the shifting spectrum of the 
source coming in and out of the detector bandpass, which might 
further complicate the structure of the detected lightcurves. 

To conclude, we remark on a few potential observational im¬ 
plications. First, we note that for all reasonable values of param¬ 
eters and for all accretion rates the jet power (as long as the jet 
exists) is at least one order of magnitude larger than the disk lu¬ 
minosity (see, e.g.. Fig. 0. This is consistent with observations 
that jet power exceeds accretion disk luminosity in blazars and that 
the accretio n flow in these systems is in the magnetically-arrested 
dis k regime dRaw lings & Saundersll99ll:IZa maninasab et alj2014l : 


iGhisellini et al. 120141 : Zdziarski et alj 2015 ). Moreover, if the j 


; jet 


is relativistic its radiation would be beamed and enhanced further. 
Secondly, we note again that both the jet power (and its correspond¬ 
ing X-ray emission) and the disk luminosity do not follow in a sim¬ 
ple manner the mass accretion rate. This should be taken into ac¬ 
count when searching observationally for TDEs or when analysing 
the observed light curves_o_f TDE can didat es. Particularly interest¬ 
ing is the possibility (Krolik & Piran 2012) that some of the above 
mentioned transitions and in particular the transition from super- 
to sub- Eddington accretion or the formation of a thin disk could 
be identified. This would provide significant new independent in¬ 
formation on the parameters of the TDE and in particular on the 
masses of massive black hole. 
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